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Abstract 

We investigate the asymptotic behavior (1 < p < oo) of a per- 
turbation of a Lax or overcompressive type shock wave solution to a 
system of conservation law in one dimension. The system of the equa- 
tions can be strictly parabolic, or have real viscosity matrix (partially 
parabolic, e.g., compressible Navier-Stokes equations or equations of 
Magnetohydrodynamics). We use known pointwise Green function 
bounds for the linearized equation around the shock to show that the 
perturbation of such a solution can be decomposed into a part corre- 
sponding to shift in shock position or shape, a part which is the sum 
of diffusion waves, i.e., the solutions to a viscous Burger's equation, 
conserving the initial mass and convecting away from the shock profile 
in outgoing modes, and another part which is more rapidly decaying 
in LP. 

1 Introduction 

Consider the system of conservation laws with viscosity: 
(1.1) ut + f{u)x = iy{B{u)ux)x 

with u G is the conserved quantity, and z/ is a constant measuring trans- 
port effects (e.g. viscosity or heat conduction). As we are not considering 
the vanishing- viscosity limit — > 0, we can assume = 1. An important 
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class of the solutions for (jLlj) are the viscous shock wave solutions, i.e., so- 
lutions in the form u{x,t) = u{x — st), where the constant s is the velocity 
of the shock, and where u connects the endstates u± = 'S(iboo). With a 
simple change of coordinates, we can assume that s = (a stationary shock 
solution), u is assumed to be an element of a smooth manifold {w'^j^giRd, 
which consists of stationary solutions of (jl.lj) connecting the same endstates 
U- and u+, and u = vP. Loosely stated, we prove that, with u a solution of 
(jl.lj) and a small perturbation of n, there is a small 60 such that 

(1.2) n(x, t) - u^'> (x) = vix, t) + ifix, t) + -Q^S{t) 

where 

— -(1—-)—- 

1. !?;(•, t)|LP ~ (l+t) ^ p . By choosing appropriate (5o and the mass 
carried hy we also obtain zero initial mass for v, i.e., 
f^^v{x,0)dx = 0. 

2. is a summation of convecting diffusion waves, i.e., self similar so- 
lutions to the viscous Burger's equation with appropriate coefficients, 

propagating away from the shock and preserving the initial mass in 

--(1--) 

the outgoing modes, \ip{-,t)\LP ~ (1 + t) ^ ^ , and 

3. \5{t)\ ~ (1 -|- t)~2~^'^ and 6{0) = 0. One can view 6{t) as indexing the 
"instantaneous" shock location and shape: employing Taylor's expan- 
sion gives us 

(1.3) - = -T^^t) + 0(|(5(t)pe-'=l^l) 

which shows that ^^5{t) corresponds to a shift in the shock location 
and or shape (up to an error of order |5(t)pe~'^l^l, which decays faster 
than V in any norm). 

For the exact definitions and conditions, see the subsequent sections; espe- 
cially see Theorem 14.141 corollaries 14.151 14.161 and their counterparts in the 
real viscosity case, which comprise the main results of this paper. 

To prove the above statements, we use p.2p and initial equations for u 
and u to obtain 

(1.4) vt- Lv = TZ{v,Vx.,-Q^S{t),(p) 
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where L stands for the hnearized operator around u and IZ \s a. remain- 
der we get applying Taylor's expansion. If G{x,t;y) is the Green function 
corresponding to dt — L, then applying Duhamel's principle yields: 



we try, then, to use a continuous induction to prove the desired rates of 
decay for v. 

The observation that a perturbation of a shock wave solution to can 
be decomposed into a sum of diffusion waves and a more rapidly decaying 
term is due to T.P. Liu (see [ Liull ILiu2]). He proved the result for weak 
shocks and with the viscosity matrix B being identity matrix. To prove 
the result, he first constructed an approximate Green function using heavily 
the weak shock wave assumption and the identity matrix, and then used an 
elaborate pointwise nonlinear iteration scheme. 

We, on the other hand, have already at our disposal the Green func- 
tion bounds we need. These sharp bounds are the result of a "dynamical 
system" approach based on Evans function and inverse Laplace transform 
techniques. This approach began for the strictly parabolic case (|^Z, ZH]) 
and then was extended to many other, more physical, regimes, such as real 
viscosity case ( [.MaZ.lj - [MaZ.4| : see also [.Zl., .Z.2j ). In these papers the 
LP, 1 < p < oo, asymptotic stability and asymptotic boundedness, of Lax 
type shock profiles were stated and proved by finding sharp Green function 
bounds for the linearized equation; some hints have also been given about 
Green function bounds for the overcompressive case. This approach does not 
require any assumption of the weakness of the shock profile, and the struc- 
tural and technical assumptions made about the equation and the wave are 
rather minimal. However, less information than what Liu's approach yields 
has been obtained about the behavior of the perturbation. 

Using the very same Green distribution bounds, we prove the results 
which Liu first observed with fewer assumptions: the shock profile can have 
small or large amplitude, be of Lax or overcompressive type (to our knowl- 
edge, this result is the first rather complete result about asymptotic behavior 
of a perturbed overcompressive shock); the viscosity term {B{u)ux)x can be 
strictly hyperbolic (section^ or have the block structure of the real viscosity 
case (section [SJ . Also no pointwise bounds on initial data are required, only 
bounds on norm and moments. In return for localization of the initial 



(1.5) 





G{x, t - s; y)TZ{v, -^6{t), ip)[y, s)dyds 
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data, Liu obtains pointwise bounds for the perturbation. We only assume 
smallness of the initial perturbation and its moment in some spaces, but 
then no pointwise information is obtained (we believe, however, that with 
a similar method and some more work, and with localization of the initial 
data, we can achieve pointwise bounds similar to those Liu has obtained). 

Plan of the paper: In computing the bounds for v using ()1.5|) . we 
frequently use Young-Hausdorf inequality and norms of the different 
components of IZ. However, a term that does not yield the necessary es- 
timates this way occurs and that is when part of the Green function that 
is like a convecting heat kernel is convoluted against (ipf)y{y, s), with ipi 
a diffusion wave convecting at a different speed. Sharp estimation of such 
terms was first treated by Liu |Liulj . Here, we extract the essential fea- 
tures of his argument, to establish that similar bounds hold whenever the 
derivative of G along characteristic directions decays more rapidly than Gy : 
in particular, for the more general Green function terms we consider here. 
Most of section |2 is to compute pointwise bounds for this part of the cal- 
culations. In section El we consider the case, already well established (see 
|LZe|, IKawH IKaw2|, ICLj ). of a strictly parabolic system with a pertur- 
bation of a constant state solution. The calculations foreshadow the more 
difficult case of shock waves. In sections H] and we consider the pertur- 
bation of a shock wave solution in the strictly parabolic case and the real 
viscosity case, respectively, which are very similar; the main difference is 
that while in section ^ we use strict parabolicity to establish short time 
estimates and thereby find good bounds for Vx, in the real viscosity case 
we have to go through a long list of energy estimates to find the bounds 
we need for the derivatives. These bounds generalize similar energy esti- 
mates obtained in |MaZ.2l IMaZ.4l IZ.21 \Z.3] . which in turn generalize the 
important estimates obtained by Kawashima and others (see )Kaw | and 
references therein) for perturbations of constant states. 

2 Some preliminary computations 

In this section, we establish some pointwise bounds for the solution u{x, t) 



Here and elsewhere in this article K{x,t) and g{x,t) both denote the heat 



of 





for t > 



1 —X' 

kernel: g{x,t) = K{x,t) = {ATTt)~2e~ 
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Using Duhamel's principle we obtain from ()2.ip . 



u{x,t) = I / g{x -y,t- s){K{y - s,sf)ydyds 



(2.2) 



oo 

/•t f+ca 

/ / 9yix-y,t- s)K{y- s,sfdyds 

Jo J~oo 



The following bounds for u are essential for obtaining bounds in subse- 
quent sections. They are similar in nature to the bounds given by [Liulj . 
but the proof given here is different from that of Liu, and is somewhat more 
general. 

Proposition 2.1. Let u{x,t) be the solution of \2. 1\) given by i2.^) . If 
t > 1, then we have 

\u{x, t)\ < Ct-l {g{x, 4t) + g{x - t, At)) 
+CX{^t<x<t-Vt}(t'^^''^ +t-Ht-x)~^). 

where x stands for the indicator function, and C is a constant independent 
of t and X. 

The same result holds if in \2.1\) and \2.^) is replaced with Kx ■ 

Remark 2.2. The result just mentioned is achieved by detecting crucial 
cancelation in the calculations. In fact, if we replace the integrands in 1)2. 2|) 
by their absolute value, we will obtain the following bounds instead (See 
IHZl): 

\u{x,t)\ <C{g{xAt)+g{x-t,At)) +Cx|^<^<i_^|(a;-^(t-x)-^). 

Proof of Provosition \2.1\ As K'^ ~ K^, the proof will be stated only for . 
It would be straightforward to observe that the same argument works for 
Kx at every step. 

We begin the proof by first stating a simple lemma: 

Lemma 2.3. // < s < yi, then e « < Ce st with C independent 
of t, s and x. 

Proof. The statement of the lemma is equivalent to 



4i - 8t 

for some D, which (after some calculation) in its turn is equivalent to 
(x ± 2s)2 - 2s2 > -SDt, which holds for L» > i since s"^ <t. □ 
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In the proof of the proposition, we will, in several places, use the following 
lemma, which is due to P. Howard | Ho |. 

Lemma 2.4 (P. Howard). Let f{cr) > be a nonincreasing function on 
and /(O) < oo. Assume further that there exist constants 7 > 0,Ci; > 1 
so i/iai /(cr) > 7e^t(^~;j) on]R+. Then for a, z > 0, 



Remark 2.5. This result is sharp down to scale a" 2. In practice, this is 
often augmented with bounds obtained by other means. See for example 
lemma 5 of [HZJ or remark [4. 171 of this paper. 

In what follows we will frequently use the following properties of the heat 
kernel g, which are easy to prove: 




Proof. See page 102 of [Ho| . 



□ 



(2.4) 




(2.5) 
(2.6) 
(2.7) 



9x{x,t)\ < Ct^2g{x,2t) 

\gt{x,t)\ < Ct~^g{x,2t) 
\g{x,t)\ < Ct--2. 



Rewriting ()2.2() . we have: 



(2.8) 




) )ydyds 
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(/) and (///) are easy to estimate: 

'"v^ r+oo 



(2.9) 
(2.10) 



J-oo 



10 J-oo 

By (|^ and the above is, 

"Vt f + OO 



g{x -y,t- s){g{y - s, s) )y dy ds 



9y{x -y,t- s)g{y - s.sfdyds 



r f -1 -1 

<C / (t-s) 2s 2g{x -y,2{t- s))g{y - s,2s)dyds 

Jo J-oo 



'0 J-oo 

which is, by 



<C {t - sy2s-2g[x - s,2t)ds. 
Jo 

Now, using lemma 1231 the above is 

f^' 1 1 

< Cg{x,4:t) / {t - s)~2s-2 ds 

Jo 

/■^ 1 

< Ct 2g[x,At) / S 2 ds 

Jo 

< Crig{x,'it). 

Part (///) in (|2.8|) can be handled similarly. 
The more difficult part is part (/I) of 1)2. 8() : 



(2.11) 



II 



t-Vt r+ 



Vt 



g{x -y,t- s){g{y - s, sf)y dy ds. 



In order to estimate (//), let us write g = g{x,T). We have then, 

g{x -y,t- s){g{y - s,sf )y = {g{x -y,t- s)g{y - s,sf)s 
(2.12) - gr{x - y,t - s)g{y - s,sf 

+ 9{x-y,t- s){g\{y - s,s). 



We will do the estimates piece by piece. 
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The first part of (|2.12|) can be estimated as follows 

ft~\ft /■+00 

(2.13) 



(2.14) 



(2.15) 



(5(2; -y,t- s)g{y - s, sf)s dy ds 

/+00 
g{x - y,Vt)g{y ~-t + Vt,t - Vtf dy 
-00 

/+00 
g{x -y,t-Vi) g{y - Vi, Vtf dy. 
-00 



Using (|n)) and (fT7)) . we will have: 

"+00 



/+00 ^ 
g{x -y,Vt)g{y-t + Vi,t- Vtf dy < Cr^g^x -t + Vi, t), 

and 

/+00 ^ 
gix -y,t-Vi) g{y - Vi, Vtf dy < Crig^x - Vt, t), 
-00 

but by lemma 1221 

g{x - Vi,t) < g{x,2t) 

g{x-t + Vt, t) < g{x - t, 2t). 

These terms fit in the right hand side of ()2.3|) . 

For the other parts in (|2.12jl . we use (|2.6|) to obtain: 

rt—\/i c+00 

9t{x -y,t- s)g{y - s, sf dy ds 

\Jt J —00 
t-Vt 



(2.16) 



/■t-vt ^ 

</ s~2(t - s)-ic,(x - s,2i)<is 

J^rt 

{t — s)~^g{x — s, 2t)ds 



Vt 



+ / s^^t - sy'^gix - s,2t)ds 



=:A + B. 
If < \/t, then 



(2.17) 



A<Ct~^g{x-Vt,2t) I ds 

Jvt 

< Cr^ g{x - Vt, 2t) 

< Ct~^g{x,4:t). 
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Similarly when x > t — \ft^ we obtain A<Ct ^g{x — t, At). 



Now for \/? < X 


< 1 we have: 


(2.18) 


A = / s-2{t 


(2.19) 


<r^ 


(2.20) 


=ct-- n 

JVt 


(2.21) 


Vt 



2 /S 
t 



2 ds 

g 8 W t-" 



2 ~i2 

cr 2e da 



„ 3 ,X, 1 



Now we use Howard's lemma, lemma ITU to find that 1)2.21(1 is indeed, 
(2.22) 

(2.23) =Cr^x-^. 

If, on the other hand, | < x < t — \/t, then clearly A is majorized by 
the value already computed for x = |, of 



(2.24) 



A < Cr^ < Cr^x-^. 



also acceptable. 

Part B in (|2.16|1 can be estimated similarly. 
Now remains the last part of ()2.12|) . i.e.. 



t-Vi r+ 



Vt 



g{x -y,t- s){g^)r{y - s, s) dy ds 



which can easily be shown to be 

r-t-Vt 



ft-Wt g 

(2.25) < / s~2g{x - s,2t)ds. 

JVt 

If X < Vt, then ^^n^ is 

rt~Vt 3 
<C g{x- Vi,2t) / s"2 ds 
JVt 
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< Ct"3 g{x - Vi, 2t) 

For X > t-Vi we use a similar method. 

For Vi < X < t - we use a similar method to what we used for 
the previous case, and again invoke the result in lemma ITU to conclude that 
(|T^ < x-i 

This completes our proof. □ 

As an straightforward consequence to the above proposition we have: 

Corollary 2.6. For u{x, t) the solution of h2.1\) given by and for 

t > 1, we have: 



where \ • |lp stands for the norm in LP(M), 1 < p < +oo, and C is a constant. 
The same result holds if in \2.1\) and \2.2^) is replaced with Kx ■ 

Bounds obtained in corollary ()2.6|) will be used when working to prove 
bounds in the system case. In that setting, however, we will usually have 
a convecting diffusion wave instead of K{x, t) in the source, and part of a 
Green function, convecting at a different speed, in place of g{x,t) in (|2.2j) . 
The following result deals with those cases. 

Corollary 2.7. In \2.^) . if we replace K{y— at, t) with (f){y , t) , and g{x—y,t) 
with G{x, t; y) in h2.2\) . we will have similar hounds for u{x, t) obtained in 
proposition \2.1\) and corollary \2. provided (j) and G satisfy the following 
hounds: 



(2.26) 



(2.27) 
(2.28) 
(2.29) 
(2.30) 

(2.31) 
(2.32) 



Gy{x,t-y)\ <Gt-^g{x-y-at,2l3t) 



Gt{x,t;y)\ < Ct~^g{x-y-at,2/3t) 



\(l){x,t)\<Gg{x-bt,(3t) 
\4>y{x,t)\ <Ct-^g{x -ht,2j3t) 
\<Pt{x,t;y)\ < Ct~^g{x-bt,2pt) 



\G{x,t;y)\ <Gg{x-y-at,f3t) 



for some a^b, and some constants C, /3 > 0. 
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Proof. A review of the proofs just carried out will show that the above 
bounds are the only ones used in the proof. Hence everything works in the 
same way as before. □ 

Example 2.8. 

/- 1 -a:^ 



solves 
(2.33) 



k - P(l)xx = -{r)x for t>0 
'){x,0)=m5o t = 0, 

where 6o stands for the Dirac distribution, and m = (j){x, t) dx. 
The function (/> satisfies the inequalities in 1)2. 7jl (see | Liu2| ). Also one can 
easily see that if one puts (j)x in place of (j)^ in above argument, then again 
one will obtain similar results, as (j)x ~ 4^'^- This function <f) will be the 
prototype of the diffusion waves, which we are going to define and use in 
the next sections. 

Now that we are in the mood of working with the heat kernels, let us 
state some lemmas that we will need when dealing with systems. The proofs 
are easy and left to the reader. 

Lemma 2.9. // oi / 02 and /3i,/32 > 0, then \K{x — ait, Pit)K{x — 
a2t,P2t)\Lp < Ce^^*, for some tj > 0. 

The same result holds if one replaces K with (p from example 

And the following lemmas, which will be needed for shock wave cases: 

Lemma 2.10. Assume a > (respectively, a < 0), h{x) is a bounded func- 
tion and h{x) = 0(e~l^'l) as x ^ +00 (respectively, as x ^ —00). Then 
\h{x)K{x — at,t)\LP = 0(e~''*) for some ?? > 0. 

Lemma 2.11. Assume a,h are both of the same sign, then 
ft f+00 



/■t /■ + CCJ 

I I g{y + a{t- s),t- s)g{y-bs,sfdyds = 0{e~ 
Jo J-00 

for some r] > 0; // a and b are of the different sign then, 

[■t r+co 
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g{y + a{t- s),t- s)g{y - bs, s) dy ds = 0{t 
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Lemma 2.12. For t > and for a / 0, 

/+00 
g{y + at, t)e^-\y\Uy = 0(^-3 e"''*). 
"OO 

3 System of conservation laws with constant back- 
ground solution 

Now consider the system of conservation laws 

(3.1) ut + /(n)a. = {B{u)u^)^ 

with the solution u = u{x,t) £ M", a perturbation of the constant back- 
ground solution u = uq. 

Let u = u — u, and use Taylor's expansion to obtain 

(3.2) Ut + Aux - Buxx = -{T{u,u))x + 0{\u\^)x + 0{\u\\ux\)x 
where A = df{u), B = B{u) and T = \d'^f{u). 

Some basic assumptions have to be made: we assume f,B £ C^, df{u) 
is strictly hyperbolic, Rea{B) > and finally stability criterion of Majda 
and Pego | .Kaw,. .MP J: Rea{—ikA — k^B) < —6k^ for all real k and some 

e>o. 

Let ai, • • • , On be the eigenvalues of A (necessarily disjoint by the strict 
hyperbolicity of A), with corresponding right eigenvectors ri,--- ,r„, and 
left eigenvectors /i, • • • ,ln, normalized so that k ■ rj = 5ij. Define F*-^ and 6*- 
to be the constant coefficients satisfying 



(3.3) r(r,-,rfc) = ^r},r„ Sr, = ^ 6;. 



n 

i=l i=l 



and set (3i = b\ and 7^ = (notice that it follows from our assumptions 
about A and B that /3i > 0). Define diffusion wave in direction: ip^{x,t) 
to be the solution of 



(3.4) 



^\ + ai^x - Pi^ix = -li{^'^)x for t > -1 
(p^{x, —1) = mi6o t = —1 
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where 60 is the Dirac distribution, and rrii is the amount of mass u{x, 0) dx 
in the direction rj. See Example H2.8() . Assuming m < Eq we wih have: 



(3.5) 



9" 



dx 



LP 



i.e., has bounds hke a heat kernel. It is not difficult to observe that ip 
acts like a convecting heat kernel; especially for our interest is the fact that 
it satisfies the bounds Finally set 



Let V := u — if, hence 
(3.6) u = u + if + V. 

Notice that 



(3.7) / Vodx = 0. 



hoo 



00 



Set Vo{x) := f^^ vodx. 

Lemma 3.1. |Vo|li = {"^oldx < \xv{x,0)\i^i , assuming that the latter 
quantity is bounded. 

Proof. We can assume, without loss of generality, that vq is scalar, i.e., 
vo{x) G M. If xo < 0, then \Vo{xo)\ = \ fy^vo{x)dx\ < ^ \x\\vo{x)\dx, 
hence |xoVb(x)| < f^'^\x\\vo{x)\dx, which approaches 0, as xq —00. 
Likewise for xq > 0, we have Vb(xo) = J^'^vo{x)dx = — f^^°° vo{x)dx (be- 
cause of ()3.7p ). and a similar argument shows that |3;oVb(x)| approaches 
as Xo +00. Now assume that (ai, 02) is an interval on which Vq does not 
change sign (suppose, without loss of generality, it is positive on this inter- 
val), and Vb(ai) = and Vb(a2) = 0. Then aiVb(ai) = and a2Vo{a2) = 
(in the case ai or 02 is ±00, the aforementioned argument would work). In 
this interval, then, we will have 

|Vo(x)(ix| = / Vo{x)dx 

which, by integration by parts, is equal to 

/■02 roi2 
02^0(^2) — oiVb(ai) — / xvo{x)dx < / \xvo{x)\dx. 

J ai J ai 

Now take summation over all the intervals in the aforesaid form. □ 
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Substituting u with u + 99 in H3.2|) . we get: 



Vt + Avx - BVxx 

= -{ipt + A(fx- B(fxx + T{(f,(f)x) 
(0.0) „ , 

+ 0{\v\^ + \ip\\v\ + \{ip + v){ip + v)x\ + \ip + vf). 



Using and we get 

(3.9) 



XX ' « 



(The whole point is that, this way, we get rid of the terms (y?*)^^^ and (Pxx'''i-) 
The Green function for the hnear part of (|3.8|) . i.e., for vt + Avx — Bvxx, 

is 

G{x, t- y) = Y,i^^t)-h^^^i^rill + R{x, t; y), 
1=1 

with the remainder R, 



Rix,t;y) = Oiil + tr^^e- 



i=l 

(for proof see |LZe| ). 

Using Duhamel's principle, 

/+00 
G{x,t;y)vo{y)dy 
-00 

rt f+00 

+ / / G{x,t- s;y){J^{v,ip)y{y,s) + ^{y,s))dyds 

Jo J ~oo 

/+00 
Gy{x,t-y)Vo{y)dy 
-00 

t f+00 




Gy{x, t - s; y)T{v, ip){y, s)dy ds 

ft f+00 

+ G{x,t- s;y)'^{y,s)ds. 

Jo J ~oo 
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Theorem 3.2. Assume the above setting, u = v + ip and assume that 
I^oIlI) |'Uo|l°°5 \xuq{x)\ii < Eq, for sufficiently small Eq (these inequalities 
translate into similar ones for vq and ipo). Then, 

(3.11) \v{;t)\p<CEo{l + ty'^^^-p^-'^ 
for some constant C. 

Proof. Let 

(3.12) C{t):= sup \v{;s)\lp{1 + s)-^^^--^^+K 

0<s<t, l<p<oo 

i.e., \v{-,s)\lp < (1 + s) ^ ^ '*C(*)) ^iid in particular \v{-,s)\l'>° < (1 + 

3 

s)~i({t). The goal is to show: 

at)<c{Eo+atf). 

But then, if Eq is sufficiently small, this implies that C(t) < 2CEq, and that 
is what we are looking for. 

Obviously \Gy\LP < r^^^"p^"i Also \ip{-,t)\LP < ^o(l + and 

— 1(1— I)— I 

\^x{-,t)\LP < -£^0(1 +0 ^ p' 2 _ need to find some bounds for Vx, and 
that is the subject of the following lemma: 



Lemma 3.3. Given above setting, we will have: 
(3.13) \vxi-,t)\LP< 



Ci\vi-,t-l)\LP + \Miipi;t)\Lp), fort>l 
C{t-^vo\Lp\ + \M{ip{-,t)\Lp), fort<l. 



where M{(p) ■= — ft + (-^(ti + (p)fx)x — {f{u + f))x o.nd consequently 

\c(B„ + ((t))t-i fort<l. 

Proof. 

ut + f{u)x = {B{u)Ux)x 

implies 

Vt + ^Pt + f{u + (p + V)x = {B{u + (p + v){Ux + fx + Vx))x- 

Therefore, 

vt + {f{u + (p + v) - f{u + ip))x - {{B{u + (p + v) - B{u + (p)){ux + y:>x))x 
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-{B{u + ip + v)Vx)x 
= -(ft - f{u + ip)x + {B{u + ip)ipx)x 

=:M{v). 

From here we can use short time estimates described in [ZHJ, section 11, 
to achieve inequahty (|3.13|) (see the argument in lemma |4.13|) . (|3.14j) fol- 
lows immediately, if we notice that, by definition of ip, we have M{'-p) = 

o{\v\Wx\). □ 



Returning to the proof of the theorem, whenever < s < t, then. 



(3.15) 



\Hy: ^){; s)\l. < c{Eo + e{t))s-"^^'-"vy 



(we used Eq < 1 and so E^,E'^ < Eq and EoCit) < \eI + \C{t) ). 
When f < 1, then 



H-,t)\Lv < C\vo 



LP 



+ 



\Gy\ii\J^{v,ip)\Lpis)ds 



(3.16) 



+ / {\G\Limy,s)\Lpis)ds 







<CEo + iEo + Citf) I (t- s)-^2s-^2ds + CE, 
Jo 

< ciEo + citf) < c{Eo + c{tf)ii + ty^- 

For t > 1 we use again Haussdorf- Young inequality to obtain: 

/+00 



'0 [\l + s)-^ 
Jo 



Gy{x,t;y)Vo{y)dy 



LP 



< I^oIliICJlp <C^ot"5(i-i/p)-i 
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and 



t /■ + 00 



J-oo 
i/2 



(3.18) 



Gy{x,t- s;y)J='{v,(f){y,s)dy ds\LP 

rt/Z 

< / \Gy\Lp\J^{v,ip){-,s)\Lids 
Jo 

+ / \Gy\LA^iv^'^)i-is)\Lpds 

Jt/2 



t/2 



Jt/2 



ids) 



It remains only to deal with the term ^ in ()3.8() . which, by (|3.9j) . includes 

the terms in the form T'^-j^{(^ ip^)xri for j 7^ A;, and r^jj{if>^)'^ri and bljipixri 
for i / J. 

Lemma [2. 91 takes care of the terms in the form {ip^if^) for j 7^ /c, as strict 
hyperbolicity of A implies aj 7^ flfc., hence giving us 



(3.19) 



1/1 1 \ 3 



This then will be treated similar to the way ()3.15|) is treated in 1)3. 18(1 . 
For other terms, we need to estimate 



t r+00 



G{x,t-y){^)y 



J-oo 



LP 



for i 7^ j. As the remainder in G is small enough, the only part of concern 
would be: 

/ (47rt)-2e 4/3,i (-^j)2 
' ./ — 00 /^p 

for i ^ j. Here we have a heat kernel convecting at the speed Oj convoluted 
against a diffusion wave (/?-^ which is similar to (j) ™- Example 1)2. 8(1 . but 
convecting at the speed Oj. Corollary (|2.7|) . then, implies that the above 

term would be less than CEq{1 + t) '^^ p' ^ . 

The term bljipxxVi can be treated similarly. This finishes the proof. □ 
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Remark 3.4. In the scalar case, we get rid of the terms in the form of 
ip^ all together, since z = j = 1 is the only possibility. Hence it can readily 
be seen that, in the scalar case, we would obtain, using this argument, the 
decay rates: \v{-,t)\LP ~ (1 + t) p' 2 ^ for e arbitrarily small | Liul |. 

4 Strictly parabolic cases with a viscous shock so- 
lution 

We now focus on a shock wave solution of the system of viscous conservation 
laws 

(4.1) ut + f{u)^ = {B{u)u^)^, 

where / € M" and B{u) G M"^", and u G M" is a perturbation of (without 
loss of generality) a stationary viscous shock solution 

(4.2) u = u{x), lim u{x) =: u±, 
i.e., u solves 

(4.3) u' = Biur\f{u)-f{u^)). 

Following |ZHl IZ1 |. we make assumptions (Ti) below. 
Assumptions (Ti.): 

(m) f,BeC^. 

(ni) Rea{B)>0. 

{Ti.2) cr{df{u±)) real, distinct, and nonzero. 

(H3) Rea{-ikdf{u±) - k'^B{u±)) < -9k^ for all real k, some 61 > 0. 

(TC'i) All set of the stationary solutions near u of (|4.1|) - (|4.2|) . connecting 
the same values u± forms a smooth manifold {u^},6 G M.^,u^ = u. Moreover 
the stable manifold of ii_ and the unstable manifold of ti+ (with respect to 
1)4. 3|) ') are transverse. 

Condition (7^3) is the stable viscosity matrix criterion of Majda and 
Pego, corresponding to linearized stability of the constant solutions u = u± 
|MP. Kawj (clearly necessary for stability of u{-) of the type we seek, see 
further discussion f |ZH |. pp. 746, 767, and 774-775). Note that condition 
(7Y4) is the condition HA of [ZHJ, plus the assertion that the shock is of 
"standard" or "pure" type (see [ZHJ, section 10). This implies that we 
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have n + i incoming characteristics, entering the shock, hence n — i outgoing 
modes, i.e., eigenvalues are in the form: 

(4.4) Oj^ < • • • < OpL^ < < Op < • • • < a„ , 
and 

(4.5) at <■■■ < a++^_i < < a+_^^ < • • • < a+, 

where af denote the (ordered) eigenvalues of df{u±). If ^ = 1, we have a Lax 
type shock wave, in which case there are n—1 outgoing modes (corresponding 
to af ^ 0) and n + 1 incoming modes (corresponding to ^ 0). If ^ > 1, 
we have an overcompressive shock, with n — i outgoing modes and n + i 
incoming modes. For further discussion see jZH| . 

The following Lemma proved in [MP| asserts that u± are hyperbolic also 
in the ODE sense (for an alternative proof, see Remark 2.3 in section 2). 
This implies exponential approach of to its asymptotic states at x = ±00, 
a fact that will be crucial in our subsequent analysis. See |MP| | and also 
[ZHJ for proofs. 

Lemma 4.1. Given (TCO) — (TCS), the stable/unstable manifolds of df{u±) 
and B{u±)^^df{u±) have equal dimensions. In particular, B{u±)~^df{u±) 
has no center manifold. 

Corollary 4.2. Given {HO) - (HA), solutions of ||]^ are in C^, satis- 
fying 

DiD\{u\x) - u±) = 0(e-"l^l), a > 0, < j < 4, i = 0, 1, 
as X — > ±00 

Linearizing about n(-) gives: 

(4.6) vt = Lv := -{Av)^ + {Bv,^)^, 
with 

(4.7) B{x) := B{u{x)), A{x)v := df{u{x))v - dB{u{x))vu^. 

Denoting := A{+oo), B^ := B{+oo), and considering corollarv 14.21 it 
follows that 

(4.8) |^(x) -^"1 = 0(e-''l^l), |B(a;) = O(e-^l^l) 
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as X — > — oo, for some positive r]. Similarly for A'^ and , as x — > +C!0. 
Also \A{x) - A^\ and \B{x) - are bounded for all x. 

Define the (scalar) characteristic speeds af < ■ ■ ■ < a^ (as above) to be 
the eigenvalues of A'^, and the left and right (scalar) characteristic modes l^, 



+ 



to be corresponding left and right eigenvectors, respectively (i.e., A^i 



afff, etc.), normalized so that l~j • r'^ = 51. and I- ■ r^^ =5;^. Following 
Kawashima |Kaw | , define associated effective scalar diffusion rates : j = 
1, • • • , n by relation 



(4.9) 



diag L^B^R^, 



\o Pi J 

where L"^ := {if, . . . ,Z^)*, := (rf , . . . ,r^) diagonalize A"^. 
Let 

(4.10) G{x,t;y) ■.= e^'5y{x) 

be the Green's function associated with (dt — L). Then, the relevant lin- 
earized theory can be summarized in the following two propositions, proved 
in EHl. 



Proposition 4.3. Given (7i ), necessary conditions for -linearized orbital 
stability, p > 0, of u{ ) with respect to perturbations G Cq" are: 
Assumptions (V): 

(PI) L has no (Lp' , without loss of generality) eigenvalues in {ReX > 
0}\{0}. 

{V2) {r±; a± ^ 0} U {/+^ ^dx; i = l,--- ,i} is a basis for W, with 
^^dx computed at 6 = 0. 

Proposition 4.4. Under assumptions (TC), (T>), we have for y < the 
decomposition 

(4.11) G = E + S + R, 

where 

E[x,t-y) ■.=j:Y. \<^:Mi-K' ( ^ 



errfn 
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(4.13) 

S{x,t;y) ■.= X{t>i} E r,-/,-*(47r/3,-0-^/'e-(--^-»^"*)VW* 

a^<0 
«fc >o 



Oj. >o, a. <o 



ifc >0,aj>0 



gX _|_ g X 



+ e 



with 



(4.14) 4(y,t):=a±(t 



and 



(4.15) ^ji(.,,,):=|.,^ + J|LgV,-, 



and 
(4.16) 
R{x,t;y) = 

+ ^ O ((t + l)-V2e-'?^+ + g-^l^l) ^-l/2^-(x-y-a^tf/Mt 
k 

+ E X{|a,-*|>|.|}0((i + l)-^/^t-^/^e-(— 7(*-l-K-|))V^*e-'^^^ 

a;^>0,aj<0 

-l2 



+ E X{|a,-.|>|.|}0((t + l)-^/W2e-(— T(*-I^K-|))VMV 



a;^>0,at>0 



for some rj, M > Q, where denotes the positive/negative part of x, indi- 
cator function x^\a~t\>\y\} onefor\a'j^t\ > \y\ and zero otherwise, indicator 
function X{t>i} one for t > 1 and zero otherwise, and scattering coeffi- 
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dents [c°''_]; [cj^^] are constant, with 

(4.17) E Kr-]-7 + E i^t-n + Et^-] / S-^- 

oj<0 ot>0 »=1 



/or eac/i (note: uniquely determined, by condition (D2)), and 
(4.18) 



- - - + + r'^du^ 1 



a-5 

where Sj denotes the jth standard basis element, and with always com- 
puted at 5 = 0. Likewise, we have the derivative bounds 
(4.19) 

\Rx\ = 

+ ^ O ((i + l)-l/2i-l/2g-,x+ ^ ^.i/2^-i,-y-a-t)yMt 
k 

flfc >0,oT<0 

+ E X{|.^-,|>|,|}0((i + i)-H-v2e-(— ;(*-i^/«^i))VM*,-..-^ 

a^>0,ot>0 

(4.20) 

\Ry\ = 

+ ^ O ((i + l)-V2e-'7^+ + e-^l^l) i-^e-(^-^-"fc~*)'/^* 



+ E X{|„-.|>|,|}0((t+l)-4-V2e-(— 

flfc >0,aT<0 

+ E X{|a,-*|>|.|}0((t + l)-^t-^/^e-(— r(*-bK-|))VM*,-..-. 

flfc >0,at>0 

^ symmetric decomposition holds for y > 0. Moreover, for \x — y\/t suffi- 
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ciently large, 



(4.21) \G\<Ce — —. 

Remark 4.5. Though it was not remarked in |MaZ.3j. the terms E and 5 
are continuous at y = 0, a consequence of the respective scattering relations 
(|4.18|) and (|4.17j) . (Note that values at y = correspond to time-asymptotic 
states described by the scattering relations, which depend only on mass and 
not position of data.) 

Remark 4.6. The term e"'''^!^'"^'^*-* in R and its derivatives corrects a 
minor omission in |Z1] . This term comes from the fact that, in the far field, 
E and S decay at this rate while entire G decays at faster Gaussian rate. 
The Gaussian decay (|4.21|) was proved but not stated in [MaZ.Sj. The 
bound for is given here only for the sake of completeness, and is not 
going to play a role in our calculations. 

Remark 4.7. In [ZlJ and | MaZ.3| the above bounds have been explicitly 
stated and proved for Lax case, and only some hints are given as about the 
overcompressive case. The proof for the overcompressive case, however, is 
not very different and can be achieved following the same outline given for 
Lax case. 




Define e^, i = 1, • • • , £ for y < 

ei{y,t) := ^ [Cfc;!.]/fe* I errfn 

(4.22) 

— errfn 
Hence 

(4.23) E{x,t:y) = Y,—{x)e^{y,t) 

i=l * 

and symmetrically for y > 0. Define also 

(4.24) G = S + R. 

We have the following bounds for G and e^'s: 



y + a,t 
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Lemma 4.8. Under assumptions (Ti) and (V) there holds 

/ + 00 ^ 
Gi;t;y)f{y)dy\Lv < Cmin{|/Up, t-2(i-Vp)|/|^,}, 
-oo 

(4.26) 

/ + 00 
Gy{;t;y)f{y)dy\L. < Cmm{t-i/2| i-,(i-i/p)-i/2|^|^^|^ 
-CX) 

for allt>0,f€L^n LP, some C > 0. 
Lemma 4.9. The kernels e, 's satisfy 

(4.27) \ei,i;t)\Lr,K{;t)\L.<Ct-2(^-'/P\ 

(4.28) K,{;t)\Lv<Ct-'2('-'M-i/2^ 

for all t > 0. Moreover, for y < we have the pointwise bounds 

\eiyiy,t)\,\e^,{y,t)\ < C 2_, t~^e , 

a->0 

\euy{y,t)\ <C}2 e M^, 

a->0 

for M > sufficiently large (i.e., > 4:b±), and symmetrically for y >0. 
Proof See HI] and | MaZ.4 |. □ 
Let H solve (|4.1|) . and, using (D2), assume that 



/ + 00 ^ p 

u{x, 0) - u{x) = ^ iT^j'/'J + ^ mjr^ + / Q 



aj<0 a+>0 «=1 



9& 



with mi's and Cj's small enough. Using the Implicit Function Theorem, we 
can find such that 



/ + 00 
u{x,0) — n^°(x) = m'jrj + m^- 
-oo _ _ J. _ 



a-<0 a+>0 
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where each is just "shghtly" different from mj. Notice that this way we 



have no "mass" in any J direction anymore. Therefore, by replacing 



u 



with n^o and without loss of generality, we can assume (^o = and 

/ + 00 
u(x,0) — u{x) = rrijrj + nijr^ . 

a. <0 a+>0 

Remark 4.10. In Lax case shock waves, u^{x) = u{x + 5), hence = 
u'{x), and 5o can be explicitly computed: 5o = ci. 

Let u{x,t) = u{x,t) — u{x) and use Taylor's expansion around u^^^\x) 
to find 

(4.29) ut + {A{x)u)x - {B{x)u^)x = -{T{x){u, u))^ + Q{u, u^)^, 

where r(x)(n, n) = cf f{u){u,u) — cf B{u){u,u)ux and 

Q{u,Ux) = 0{\u\\ux\ + \uf). 

Denote = r(iboo), and note that we have similar statements to ()4.8|) for 
r(x) — r^. Define constant coefficients bfj and F^.^ to satisfy 



(4.30) T^irf,4) = Y,rf,,rt, B^rf = Y,b. 



i=l 1=1 

hence of course = bf-, and denote := F^-. 



Remark 4.11. As it is pointed out in [LiuTl ILiu2l ILiu3[ . 7^ is 

equivalent to genuine nonlinearity of the i^^ field, and jf^ = means that 
the i^^ field is linearly degenerate. 

We define diffusion waves along outgoing modes: for a~ < define the 
diffusion wave (pi to solve: 



(4.31) 



+ ai(pi- Pi = {V>' )x for t > -1 
(/j'(x, —1) = mj^o t = —1 



and likewise for af > 0, define (pi to be the solution of 
(4.32) 



^i+4^i- Pt^ix = -iHv^)^ for t>-i 



(^*(x, — 1) = m-i^O t = —1 



25 



and set 

n 

a~<0 a+>0 

Let V := u - if - ^5{t), where 6{t) = {6i{t), ■■■ , 6i{t)Y^ is to be defined 
later, assuming 6{0) = 0. Notice that 

/ + 00 
v{x,0)dx = 0, 
-oo 

so if Vq = v{y, 0)dy then by lemma IXTI Vn G and iVol^^i < IxvqIli- 
Replacing u with v + ip + ^6{t) in ( ^ computed at 6 = 60 = 0), 

93" 



and using the fact that satisfies the linear time independent equation 
= 0, we will have 



(4.34) vt-Lv = ^{x, t) + V, -g^S{t))x + -Ql^t)^ 

where 
(4.35) 

+ \i^ + v+ + V + -Q^^)x\ + W + V + 

and := -ipt - {A{x)lp)x + {B{x)ipx)x - (r(x)((^, (/?))^. For ^ we write 
(4.36) 

^'(x, t) =- {ipt + Aifx - Bifxx + r((^, <f)x) 

= - E '^i^" + - {B{x)if^,rr), + (r(x)((^Vr,^Vr)), 

a-<0 

- E '/'^^^ + ^^i^)^'4)x - {B{x)^xrt)x + (r(x)(^V+,^V+)). 

a+>0 

Let us write a typical term of the first summation {a^ < 0)in the following 
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form: 
(4.37) 

^\r- + (A(x)c^Vr), - {B{x)i/,r-)^ + (r(x)(v9Vr, ^V")), 

= [{A{x) - A-)^Wr - (Bix) - B-)^lr7 + (r(x) - r")(v9Vr, ^Vr)]^ 
+ iplrr + i^lA-r7) - ((^L^-r") + i{^%T~ {r7 , r7 )) . 

Now we use the definition of 99* in (|4.3ip and the definition of coefficients 
bij and Tijk in (|4.80j) to write the last part of (|4.37|) in the following form: 



r ■ 

II J 



Similar statements hold for af > with minus signs replaced with plus 
signs. 

Later we will need some estimates for Vx- Short time estimates gives 
us the necessary bounds. Let us first provide the requisite short time exis- 
tence/regularity theory for general quasilinear parabolic systems, using the 
paramatrix method of Levi | LSUl ILe j . 

Proposition 4.12. LetA{x,t), B{x,t), andC{x,t) be uniformly bounded in 
L°° and C^^''^^^^"'"^'^\x , t) , < a, < 1, taking values on a compact set, with 
Re(T{B) positive and bounded strictly away from zero. Then, for < t < 
T, T sufficiently small, there is a Green's function G{x,t;y, s) G C^'^{x,t) 
associated with the Cauchy problem for 

(4.39) vt = Cv + iAv)x + {Bvx)x, V G M", 
satisfying bounds 

(4.40) |D^G(x,t;2/,s)| < Ct-(-''+^)/2g-(x-j,)VM(i-s)^ j = 0,l, 

where C , M , T > depend only on the bounds on the coefficients and on 
the lower bound on Rea{B). 

Proof. See [ZHIIl5U| . □ 

The following lemma provides us with the short time estimates we need 
for Vx. 



27 



Lemma 4.13. Given the above setting, and assuming v remains bounded 
for all the time, we will have: 



C{\V{; t - l)\Lr + \M{<p, %S)\lv), fort>l 

as 



C{t-i\vo\Lp\ + IMi'p, ^S)\lv), for t < 1. 



(4.41) \vx{;t)\Lp< 
where 

(4.42) M{^, ^S) = -^S + 0{\^ + ^S\\^^ + {^6U + e-'\^\\^\). 
for some k > 

Proof ut + f{u)x = {B{u)ux)x and f{u)x = {B{u)ux)x implies 
(4.43) 

Vt + Vt + -Q^S + f{u + (p+ —5 + V):c - f{u)x 

= {B{u + ip+ —6 + v){ux + (Px + + '"x))x - {B{u)ux)x- 



Hence: 
(4.44) 



vt + {f{u + ip + —5 + v)-f{u + ip + —(5)), 



-{B{u+ip + —6 + v)vx 



= -Vt - - {f{u + (p+ ^S) - f{u))x 
+((5(u+^ + —5) - B{u))ux)x 
+{B{u+ip + -Q^5){Vx + (-Ql)^^)x 

It is not difficult to observe that 
(4.45) 

Mi^p, —S) = -<pt - {Aix)iip + -g^S% + iB{x){<px + -^S)x - —S 
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Using (lOHl) . (j07|) . (jllHl) and the fact that = 0, we conclude 

(j332I)- We use 



fiv + v)- f{r]) = / df{r] + 9v)dev 
Jo 

and similar equation for B in order to write (|4.44jl in the form 

(4.46) vt + {A{x, t)v)x - {B{x, t)vx)x = Miip, -^6). 

Now A and B depend on v. Momentarily assume v (hence A and B) is in 

C(0'0)+("'°/2)(x,t). 



(4.47) 



/+00 
Gx{x, t;y,t- T)v{y, t-T)dy 
-oo 

+ I I Gxix,t;y,s)M{y,s)dyds. 



By Duhamel's principle, where G is the Green's function for (|4.46|) . and 
using the Gx bounds of proposition 14.121 for divergence-form operators, we 
find that 

\vx{;t)\oo < C(T)-V2|^(.,t - T)U + C(r)V2|A^|^. 

In particular, for t > T, we obtain a uniform Holder (indeed, Lipshitz) 
bound on v(-,t) depending only on the norm of v{-,t — T). By the 
(standard) method of extension, we thus obtain uniform Holder continuity 
of V so long as \v\ remains bounded. 1)4.41(1 follows using (|4.47() and taking 
(without loss of generality) T = 1. For a more detailed discussion see [ZHJ 
(section 11). □ 
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Now we employ Duhamel's principle to get from H4.34() : 

/ + 00 
G{x,t;y)vo{y)dy 
-oo 

rt p+oo Q^5 

~^ J J Gi^^t- s;y)J^{ip,v,—6{s))y{y,s)dyds, 

ft /"+00 



(4.48) 



+ 



/ / G{x,t- s;y)^{y,s)dyds 

Jo J -oo 



+ 00 

Gy{x,t;y)Vo{y)dy 

-00 

t r+00 Q-5 




10 J-00 
ft /"+00 



Gy{x, t - s; y)J^{(p, v, -^6{s)){y, s)dy ds, 



+ 



/ / G{x,t- s;y)^{y,s)dyds 
Jo J -00 



(the last part of the above equation follows from G{x,t]y)^-{y)dy 
= m- and 5{0) = 0). Set 

/oo 
e^y{y,t)Vo{y)dy 
-00 

(4.49) + / / e.,y{y,t-s)J'{^,v,—5{s))dyds 

i/ J ~ CxD 

pt p+00 

- ei{x,t - s;y)^{y,s)dyds. 

Jo J-00 

Using and G = S + G we obtain: 

/+00 
Gy{x,t;y)VQ{y)dy 
-00 

(4.50) - j J Gy{x,t- s;y)J^{ip,v,—6){y,s)dyds, 

pt p+00 

+ G[x,t- s-y)^{y,s)dyds. 

Jo J-00 
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We are now in possession of the necessary tools to state the fohowing 
theorem: 



Theorem 4.14. Let (7i) and (D) hold, and \uo\]^i, \xuo\£^i, \uo\loo < Eq, 
Eq sufficiently small (these assumption on uq are being inherited by vq). 
Assume the above setting and u = v + Lp + then for any e, < e < i, 



(4.51) \v{-,t)\Lv <CEQ{l + t) 
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|{l-l/p)-i 



(4.52) \5{t)\<CEo{l + t)-2+', 

(4.53) \5{t)\<CEQ{l + t)~^+\ 

for any p, 1 < p < oo, and with C independent of p (but depending on e). 

An immediate consequence to this theorem is the following corollary, 
which is almost (up to an e) Liu's result. 

Corollary 4.15. 

. , , A , f(l + t)-5(i-Vp)-| forl<p<^, 

However, our approach yields more information about the behavior of the 
perturbation, as we can track the shock location: by H1.3|) and the comment 
right after, the following corollary follows. 

Corollary 4.16. 

(4.55) \u - - if\Lv < (1 + t)-5(i-i/p)-| 

for all p. 

Proof of Theorem \4.14\ Fixing e, define 
(4.56) 

C(t) := sup |u(-,s)|lp(1 + s)^(^"p^+2 + sup |<5(s)|(l + s)^-^ 

0<s<t,l<p<oo 0<s<t 



+ sup |(5(s)|(l + s 

0<s<i 



,1-e 
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Our aim is to show that 

(4.57) at)<C{Eo + e{t)) 

and then use a straightforward continuous induction. Equivalent to ()4.57l) 
is 



(4.58) \v{; s)\lp < C{Eo + C'(t))(l + s)-^(^-f) 

and 



1 \ 1 

4 



(4.59) \s^s)\<C{Eo + emi + s) 



and a similar statement for 6. We need to take three steps: 

9(5. 



step 1: bounds for |v|lp: By corollary 14. 21 we have \^-\ ~ e '^'^l and 
""xl ~ e~'^l^l for some A; > 0, hence, for < s < t we have: 



(4.60) <cc(^)(l + s)-5+^ 

3ffi 1 

(4.61) |(_),5(.)|z.. <CC(^)(l + .)-5+^ 

(4.62) |v;(-, s)|lp < C7£;o(l + sy'^^^-'i\ 

(4.63) |v5x(-,s)|lp <CSo(l + s)~^(i-i)-5. 

Lemma l4. 1 31 provides us with the necessary bounds for Vx- Note that by 
lemma [2. 101 and the above bounds, we have the following bounds for in 



(4.64) \M{^, ^S){; s)\lp < C{Eo + C(0)(1 + s)"^^'"^^"^ 



when 1 < p < 2. Therefore, 
(4.65) bx(-,s)|LP< 



C(^o + C(i))^"^^^"'^^"^> fors>l 
C{Eo + C{t))s-^, fors<l. 



Now let p* = hence ^(1 - ^) + | = 1 - 2e. With bounds for ip, 6, 
and Vx we obtain: 

(4.66) \Hv, s)\l. < CiEo + Citf)s-^(^-^M-i 
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whenever 1 <p<p*. 
When t<l, then 



\v{-,t)\LP < C\vo\lp 



f ~ du 

+ \Gy\LA^{^^'f-:-Q^^)\Lv{s)ds 

(4.67) + A|GUi||V'(2/,s)U.(s)d5 

< C£;o + (£^0 + CW^) / (t-s)-^s-^ds + C£;o / (1 + 

< c(£;o + c,{tf) < c{Eo + c(i)')(i + 

For t > 1 we use again Haussdorf- Young inequahty to obtain: 

" + 00 



(4.68) 



and 



/ + 00 
Gy{x,t; y)Vo{y)dy\LP 
-oo 



t/2 /.+00 _ Q-S 

/ Gy{x,t- s;y)J^{v,ip, —5){y,s)dyds\LP 

-/-oo 



< / |G2;|z,p|J^(i;,(/5, — 5)(-,s)|Lids 

(4.69) ° 

< ^(£^0 + C(i)^) / {t - s)~^^^-^M-h-^ 

Jo 

<<^(Co + C(i)')i"^^'"'/*'^"^ 
<2C(Co + C(i)')(l+0"^^'"'/^^"^- 



If 1 < p < p*, then 



It/2 J-oo 95 



Gy{x, t - s; y)J='{v, ip, -^S){y, s)dy ds\LP 



(4.70) 



< / \Gy\LA^{v,cp, —S){-,s)\Lpds 

pt 

<C{Eo + C{tf) / (i-s)-is-5(i-Vp)-|t 

Jt/2 

< C(£o + C(i)^)i"^(i-VpM 

< 2C{Eo + C{tf){l + t)-^^^-^/P^-'^. 
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If p> p*, then choose q so that | + 1 = + | , and then 



t/2 



J Gy{x, t - s; y)T{v, ip, -Q^S){y, s)dy ds\Lv 

r-t 



(4.71) <C{Eo + C{tf) [ (t-s)"^(^"i)"5s-5(i-^^)-!ds 

Jt/2 

< C(Co + C(i)^)i"^^^"^^^^^^ 



It remains to show 

\[ r G{x,t- s-y)-^{y,s)dyds\Lv <EQ{l + t)-'^^^-v^--i. 

Jo J -oo 

By 1)4. 36() . (14. 37(1 and ()4.38|) we have to estimate the following: First, the 
terms in the form {LpiLpjT{x){rf ,r^))x, in (|4.3()j) . By lemmaEIHl {'fi'fj'^{x){rf, 
is of order £'gO(e~''*),, so we can use similar calculations as before to con- 
clude that \ !l Gy{x,t - s;y){ipiipjT{x){r^,rf))dyds\LP = 0{Eo{l + 
t)-|(i-i/p)-i) 

Next, the terms in the second line of H4.37() . i.e., in the form ({A{x) — 
A~)(p^r~)x or similar forms, of which lemma V2.1U\ together with 1)4. 8(1 and 
similar bounds for T take care. 

Finally, the terms in the form (flxb^j^J and (¥'*)^r~..rj, i 7^ j, in (|4.38j) . 

The Lp norm of these terms is of order (l+t) p ^ , so not decaying fast 
enough to use calculations similar to what we have already done, so we have 
to use the results in section [21 Assume a~ < 0; we examine the integration 
of {^p^)'^^~ against the different terms of S{x, t — s;y). As l^^^^fj = if A; 7^ j 
the terms of concern in S integrated against ((/?*)^^rj are as following: For 
y < 0, the terms in the first line of 1)4. 13() . gives us: 

(4.72) j'^ j {AT,(i-{t- s))-^l\-^''-y-''^^'-'^^"'^^^^'-'\^,)ldyds, 

which, by and (fTTj) . is of order £^0(1 + i)~^^^~^^~^ (as i / j). The sec- 
ond line of (|4.13j) does not comprise anything. The third line of S comprises: 
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(4.73) 

with > 0, a~ < 0, and z~i^ and computed at f — s. Note in this case 
the convection and diffusion coefficients are not constant. To make a brief 
presentation of the relevant calculations, we first notice that the biggest 
part in (|4.73|) is in the cone a~{t — s)/2 < x — zjf^ < —ci~{t — s)/2, or 
equivalently, 3aJ (t — s)/2 < x — ajy/a^ < aj{t — s)/2 (outside this cone we 
have a negligible term). On this interval the diffusion coefficient can be 
bounded from above by a constant The derivatives of Pjf^{t — s) also can 
be bounded from above similarly. As a consequence the y and t derivatives 
of this part of S satisfy the bounds used in proposition 12.11 in the cone just 

mentioned. Also we make a change of coordinates z = to see that , this 
part of (|4.73() can be estimated the same way one would estimate 

g{x - z - a~{t - s),l3*{t - s))g{z - aj s/a^ , s)ldyds. 

using the same process as in proposition 12. II and subsequent results. Notice 
that a^aj /a^ ^ aj , i.e., the different speed of the Gaussian kernel in the 
Green function and the diffusion wave. 

For y > 0, the corresponding terms in second line of S in (|4.13|) gives: 

The case would be different from (|4.72|) if — ■ I'^ this case we have a 
term like: 

^'^-^^^ (^ gx^g-x ^ j J 9{x-y-a{t-s),t-s)g^{y-as,s)dyds, 
with a < 0, which, by some elementary calculations, is less than or equal to: 




(4.75) g{x-at,Mt). 

gx _|_ g X 

Now use lemma 1^.101 

The terms in the remainder R{x, t; y) can be dealt with in a similar way. 
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step 2: bounds for 5{t): In order to show H5.27() holds we investigate 
the integration of and against each term in Vo,J-{v,ip, -^5) and 
respectively. To that end, we will have the following (assume t > 1). 



+ 00 

eiy{y,t)Vo{y)dy 

-oo 

^^''^^^ <\eiy\L^\Vo\Lidy 

and 
(4.77) 

I / / eiy{y,t - s)0{\v\^){y,s)dyds\ 

Jo J-oo 

ft/2 

Jo 

+ C \eijLi{y,t - s)\\v\'^\Loo{y,s)ds 

Jt/2 

ft/2 



r'^ 1 r 3 

<CC\t){ {t-s)'Hl + srHy,s)ds+ (l + s)-2ds) 

Jo Jt/2 

<cc2(^)(l + ^)-i+^ 



Similarly for 0(1^1 For 0((^|-(5)^) we use lemma 12. 121 to get 
(4.78) 



10 J~oo 96 



eiy{y,t - s)0{\—6\ ){y,s)dyds 

ft , f+oo f)„-.S 

<Ce{t) (l + s)~2+^y e,^(y,t-5)(|^P)(y,s)d2/ds 

rt 

<CCHt) / (l + s)-^+^(t-s)-^e-^(*-*)ds 

JO 

<cc\t){i+tr-2. 

For 0(|99p) we use the fact that both (p and are the summation of signals 
like convecting heat kernels, moving away from shock. Hence using lemma 
12.111 gives us: 

/•t r+oo ^ 

(4.79) 1/ / e^y{y,t-s)0{\ip\''){y,s)dyds<CEo{l + t)-2. 

Jo J-oo 
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All the other terms in J-{v, if, ^-5) and ^' can be treated with similar meth- 
ods. 

step 3: bounds for S{t): Very similar to the previous calculations for 
6{t). □ 

Remark 4.17. As for the conditions on initial data un, it is in fact enough to 

3 

assume (as Liu and others have done) that uo = 0(1 + |2;|)~2. To see briefly 
why this works notice that the only part we should change in our argument 
is the linear part (|4.68() . Now for this part, it is enough to consider the 

3 

convolution of a Gaussian signal, say g(x,t), against vq- As |i;o| ~ (l + |x|)~2 
and I Vol ~ (1 + |x|)~2, we consider gy{x — y,t)Vo{y)dy when |x| < y^, 

and g{x — y,t)vo{y)dy when |a;| < \/t. It is not difficult to observe, 
using Howard's lemma 1231 that 

'•+«^ ~ _1 _1 _3 

G{x-y,t)vo{y)dy\ <C{x{\^^<^tyt 5(l+|x|) 5+X||^|>^|(l+|x|) 2). 

The necessary bounds then follows immediately. To obtain bounds 
in the case |x| < Vt we use lemma 1231 to see that 

+ CXD 

9y{x - y,t)Vo{y)dy 

■00 

f +00 



J —00 

Vt' Vt 



<ct-'i r^e-yy'\y ^-"^y 

J — 00 



Other bounds follow using interpolation. 

Remark 4.18. It is an easy observation that, in order to have (|4.66j) for 
1 < p < p* < 2, it is enough to have bounds for Vx in ()4.65|) . as Vx is 
always multiplied by a favorable term in This will become important 
when we consider the real viscosity case, as we have only good bounds 
(using energy estimates) for Vx in that case . 

Remark 4.19. The Analysis can go through in the case df{u±) is not 
strictly hyperbolic, provided that A± and B± are simultaneously symmetriz- 
able, by replacing Green function bounds with more general bounds given 
in proposition 5.10 of |Z.3| . and replacing the diffusion waves (/3j's with 
the "multi-mode diffusion waves" of |Ch) and jLZeJ . The same remark is 
applicable in the real viscosity case. 
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5 Real viscosity case 

In this section we follow closely the notations and assumptions used in |Z,3j . 
Consider a general system of viscous conservation laws 

Ut + F{U), = {B{U)U,),, 



(5.2) u ={:;,], B ' 



modeling flow in a compressible medium. We assume 

ui e M"-'', u^^ G and 

(5.3) Reab2 > 9, 
with 9 > 0. 

Again we consider the viscous shock wave solutions of ()5.1() . which are 
in the form: 

(5.4) U{x,t) = U{x), lim U{x) = U±, 

X— >itoo 

satisfying the traveling- wave ordinary differential equation (ODE) 

(5.5) B{U)U' = F{U) - F{U-). 
Considering the block structure of B, this can be written as: 

(5.6) F^{u^,u^^) = F^{uL,uL^) 
and 

(5.7) biiu'y + b2{u''y = F'\u',u'')-F'\uL,u'J). 

We assume that, by some invertible change of coordinates U W{U), 
possibly but not necessarily connected with a global convex entropy, followed 
if necessary by multiplication on the left by a nonsingular matrix function 
S{W), equations 1)5. 1|) may be written in the quasilinear, partially symmetric 
hyperbolic-parabolic form 



(5.8) + AW, = {BW,), + G, W 



^11 I ^ 
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g ^n-r^ ^11 ^]SJ',xe R'^, t G M, where, defining W± := W{U±): 
(Al) A{W±), := All, A^ are symmetric, i° > 0. 

(A2) No eigenvector of dF{U±) lies in the kernel of B{U±). (Equivalently, 
no eigenvector of A{A9)~'^{W±) lies in the kernel of B{W±).) 



(A3) 13 = 




, with Reh{W) > 9 for some 6* > 0, for all W, 



and g{W^, W^) = 0{\W^\^). Here, the coefficients of dEHl 

may be expressed 

in terms of the original equation (|5.1|) . the coordinate change U W{U), 
and the approximate symmetrizer S{W), as 

io := S{W)idU/dW), A := S{W)d{dU/dW), 
^^'^^ B := S{W)B{dU/dW), G = -{dSWa;)B{dU /dW)W^. 

For examples about Navier-Stokes and Magnetohydrodynamic equations, 
see |Z.2 |. Along with the above structural assumptions, we make the tech- 
nical hypotheses: 

(HO) F, B,W,Se C^ with s > 5. 

(HI) The eigenvalues of A^, are (i) distinct from 0; (ii) of common sign; and 

(iii) of constant multiplicity with respect to U . 
(H2) (T{dF{U±)) real, distinct, and nonzero. 

(H3) Local to U{-), solutions of (|5.4|1 - (|5.5|1 form a smooth manifold {0^ {•)}■, 

Analogous to lemma UTTl and corollary 14.21 we have the following lemma 
proved in [MaZ.Sj. 

Lemma 5.1. Given (H1)-(H3), the endstates U± are hyperbolic rest points 
of the ODE determined by jS. 7| ) on the r -dimensional manifold i5.6]) . i.e., 

the coefficients of the linearized equations about U^, written in local coordi- 
nates, have no center subspace. In particular, under regularity (HO), 

(5.10) ^^^^[/^(x) - C/±) = 0(e~"l^l), a>0, 0<j <6, i = 0,l, 

as X — > iboo. 

We now recall some important ideas of Kawashima et al concerning the 
smoothing effects of hyperbolic-parabolic coupling. The following results 
assert that hyperbolic effects can compensate for degenerate viscosity B, as 
depicted by the existence of a compensating matrix K. 
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Lemma 5.2. ([KSh]) Assuming , A, B symmetric, AP > 0, and B > 0, 

the genuine coupling condition (GC) No eigenvector of A lies in keri? is 

equivalent to either of: 

(Kl) There exists a smooth skew- symmetric matrix function K{A^ , A, B) 
such that 

(5.11) Re {K{A^)-^A + B){U)>^. 
(K2) For some 6 > 0, there holds 

(5.12) Re ai-iaA'^r^A - l^iA'^y^B) < -0|^|V(1 + I^P), 
for all ^ G R. 

Proof. These and other useful equivalent formulations are established in 
see also [Ol lMaZ:4llZ3| . □ 

Now returning to the original equation (|5.H) with U a shock solution we 
linearize around U exactly as we did in section |1J and we define again A and 
B in the same manner: 

(5.13) B{x) := B{U{x)), A{x)V := dF{U{x))V - dB{U{x))VU^. 
Assume for A and B the block structures: 

\A2l A22/ \i>21 -D22 

The characteristics speeds af, the left and right eigenvalues if ,rf for dF{u±) 
and f3f = ifB^rf are all defined the same way as before and again we have 
/?f >0. 

Also, let a*(x), j = 1, . . . , (n — r) denote the eigenvalues of 

^* := Ml — M2B22 

with lj{x), rj{x) € M""*" associated left and right eigenvectors, normalized 
so that Ij^rj = dj^. More generally, for an m^-fold eigenvalue, we choose 
(n — r) X m'j blocks L* and R* of eigenvectors satisfying the dynamical 
normalization 

L*^dxR*j = 0, 
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along with the usual static normalization L*^Rj = Sj^Im* ; as shown in 
Lemma 4.9, [MaZ.lJ, this may always be achieved with bounded L*, R*. 
Associated with L^, R'j, define extended, n x m* blocks 



L*\ / R* 



Eigenvalues a*j and eigenmodes C*j, TZ*j correspond, respectively, to short- 
time hyperbolic characteristic speeds and modes of propagation for the re- 
duced, hyperbolic part of degenerate system (|5.1j) . 
Define local, rrij x nij dissipation coefficients 

r?*(x) := -LfD,R*{x), j = 1, . . . , J < n - r, 

where 

D^x) : = 

M2B22 Ml — ^22-622^-621 + ^*-B22^-B21 + B22dx{B22 B21) 

is an effective dissipation analogous to the effective diffusion predicted by 
formal, Chapman-Enskog expansion in the (dual) relaxation case. 

At X = ±00, these reduce to the corresponding quantities identified by 
Zeng [Ze.ll ILZe| in her study by Fourier transform techniques of decay to 
constant solutions {u,v) = {u±,v±) of hyperbolic-parabolic systems, i.e., of 
limiting equations 

Ut = L±U :=-A±U^ + B±U,,. 

As a consequence of dissipativity, (A2), we obtain (see, e.g., | Kawl ILZel 
IMaZ.3j . or Lemma EHU 

(5.14) pf > 0, Rea{r]*^) > for ah j. 

However, note that the dynamical dissipation coefficient D^{x) does not 
agree with its static counterpart, possessing an additional term B22dx{B22 B21] 
and so we cannot conclude that 1)5.14(1 holds everywhere along the profile, 
but only at the endpoints. This is an important difference in the variable- 
coefficient case; see Remarks 1.11-1.12 of [MaZ.3| for further discussion. 

We also make the following assumptions, necessary for linear stability: 
Assumptions (D): 

{Dl) L has no (L^, without loss of generality) eigenvalues in {ReX > 
0}\{0}. 

{D2) {rf;af ^ 0} U {/+^ ^dx; i = 1, •••,£} is a basis for M", with 
I-^ ^-(ix computed at (5 = 0. 
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Proposition 5.3. \MaZ.3fj Under assumptions (A1)-(A3), (H0)-(H3), 
and (D1)-(D2), the Green distribution G{x,t;y) associated with the lin- 
earized evolution equations may be decomposed as 

G{x,t;y) = H + E + S + R, 

where, for y < 0: 

J 

H{x,t;y) := ^ a*"l(x)a*(y)7^*(x)C;(y, t)<^.-a*t(-y)/: -'(y) 
(5.15) f 

= ^7^*(x)0(e-''«*)5._,.,(-y)£f(y), 

where the averaged convection rates a* = a*{x,t) in I15.15\) denote the time- 
averages over [0,t] of a*{x) along backward characteristic paths z* = z*{x,t) 
defined by 

dz* /dt = a*{z*), z*{t)=x, 

and the dissipation matrix = C*{x,t) £ M'^j^'^i is defined by the dissipa- 
tive flow 

dq/dt = -ri*{z*)C*, C;(0) = /m,. 

E and S have exactly the same form as in proposition \4-4[ '^^^ 
(5.16) 

Rix,t;y) = 0(e-''(l^-^l+*)) 

+ ^ O ((t + l)~V2e-r,x+ ^ ^-r,\x\\ ^-l/2^-ix~y~a-t)yMt 



k=l 



E ^{|aT*|>M}0((* + l)-i/^i-V2)e-(-«7(*-I^K-|))VMt,-..-^ 



aj. >0, a. <0 



1/2.-1/2N -(x-a+(t-|j//a^ |))2/Mi -ryx" 



aj. >0,aj>0 
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(5.17) 

J 

Ry{x, t;y) = J2 0(e-''*)5._,.i (-y) + 0{e-^(\-v\+')) 
i=i 

n 

-{x-y-altf/Mt 



+ O ((t + l)-V2e-^-+ + e-^l^l) t-^e 

k=l 
a~>0, aj<0 

+ E X,.-,>|,||0((t + l)-V2t-)e-(-4(-I^KT|))VM.^ 



2;„„-f 



a^. >0, aj>0 



(5.18) 

J 

i?,(x, t; y) = 0(e~''*)5,_,*t(-y) + 0(e-^(l^-^l+*)) 



+ ^ O ((t + l)-ie-''"^ + e-^l"l j t-i(t + i)i/2e-(-y-a,-t)VMt 

k=l 

+ E X{|,-,|>|,|}0(t + l)-V2t-i)e-(-'^7(*-I^K'l))V*«e-''^-*- 



a^. >0, a. <0 



-1/2,-lN -(a;~a+(t-|s,/a^ I)) VA^*-'?!:- 



Oj. >0, aj>0 

Moreover, for \x — y\/t sufficiently large, \G\ < C'e~''*e~l^~^l^/*^*) as in the 
strictly parabolic case. 

Once again let G = S + R and define e^'s as before. Obviously tlie same 
bounds mentioned for S and e^'s in lemma hold here also. Furthermore 
we have the following for H: 

Lemma 5.4. With the conditions in vrovosition \5.!A H satisfies: 

+00 

Hi;t;y)fiy)dy\LP<Ce-^y\LP, 

-oo 

+ 00 

H,{;t;y)f{y)dy\LP<Ce-''y\wi,P, 

-oo 

for some C,r] > 0, for any p > 1 and f £ W^'^ . 
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Proof. See |MaZ.3L lMaZ:4| . □ 

From here on, almost everything would be very similar to the strictly 
parabolic case in sectional by replacing JJ with JJ^'\ for a small 60, we 
may assume that the initial perturbation has no mass at the directions, 
then we define diffusion waves, (pi, exactly as in ()4.31|) and ()4.32|) . then 
^ = ^a-<o + ^a+>o ^^ain V = U-iJ^° -if- %-d{t) with 

S to be found (from now on we once again assume, without loss of generality, 
5o = 0). The equalities H4.34() to (|4.38|) are reproduced exactly as before, 
but with u and v replaced by U and V, respectively. Furthermore, it is easy 
to see that 

(5-19) +\{V + ^ + ^^Uiv'' + ^ + 

The function 6{t) is defined as in H4.49() . and then, similar to ()4.5U|) . we have: 

f +00 



/+00 
{H + G){x,t;y)V{y,0)dy 
"OO 

(5.20) - J {H + G){x,t - s;y)J^{^,V,—6)y{y,s)dyds, 

pt f+00 

+ {H + G)[x,t- s;y)^{y,s)dyds. 

Jo J-00 

Theorem 5.5. Let (A1)-(A3) and (H0)-(H3), (D1)~(D2) hold, and 
Wo\L^r\L°°r\H^ y |xC/o|li ^ Eq, Eq sufficiently small. Assume the above set- 
ting and U = V + (f + ^§^S; then for any e, < e < |, 

(5.21) \V{;t)\Lr> <CEo{l+t)-^2(^-'/P)-k, 

(5.22) \5{t)\<CEo{l + t)-'2+', 

(5.23) \6{t)\<CEo{l + tr'+\ 

for any p, I < p < 00, and with C independent of p (but depending on e). 
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Remark 5.6. Corollaries similar to 14. 151 and 14. 16l are valid here also. 



Proof. Fixing e, define 
(5.24) 

C{t):= sup |y(-,s)|Lp(l + s)^^^"p)+^+ sup + s)5-^ 

0<s<t,l<p<oo 0<s<t 

+ sup + 

0<s<t 

To show 

(5.25) at)<CiEo + CHt)) 
we need to show 

(5.26) \V{; s)\lv < C{Eo + C\mi + s)-^(i-?)-i 
and 

(5.27) \s^s)\<C{Eo + emi + s)-h+^ 

and a similar statement for 5. From here on the proof goes very much 
similarly to the proof of theorem 14.141 except for two issues: first, here we 
do not have a lemma similar to lemma 14.131 as the short time estimates 
there need the strict parabolic hypothesis. Instead we have to use some 
energy estimates in order to control the derivatives of V. Using this method 
we will find out that, under the assumptions of the problem, 

(5.28) \V{;t)\H3 <C{Eo + C{tm + s)--^. 

This in turn will implies (|4.66j) . for 1 < p < p* , with p* as before (see 
remark f4.18() . The other difference is that, here we have the extra term H 
in the Green function decomposition, but we do not have any bounds for 
Hy. Hence we have to compute 

/ / H{x,t;y)J='y{y,s)dyds. 

Jo J~oo 

By and ((OSl), we obtain 

1 < p < p* . This with lemma EU provides us with necessary bounds. 
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It remains to show that (f??^ holds. Let IJ - U = V + if + and 

W := W-W-ip- with ip = dWy. Notice also that dW{^) = ^ 
(we use the notation dW = dW{U), dU = dU{W), etc). 

Claim: 1^1//'' ~ |1^|_H'''- 
proof of the claim: Note that U — U = dUAve{W — W), where dUAve = 
/q dU{W + e{W - W))de. Now using the facts that ^ = dtJ^ and 
if = dUip we deduce: 

- dW^ 

V = dUAve{W) + {dUAve " dU)ip + [dUAve " dU)—-5. 

00 

This, with a similar argument in the reverse direction proves our claim. 

From here on we follow closely the argument presented by Zumbrun in 
[Z.3| (see also [MaZ.4j ). with some necessary modification to handle our 
more complicated case, e.g., key cancelations in (|5.38p and (|5.39|) and the 
term ^ in (|5.48|) . which has no counterpart in fZ.3f . First, we introduce the 
weighted norms and inner product 

s 

(5.29) l/U := \a'/^f\L^, \f\Hg := l^^/U, {f,9)a ■= {af,9)L^, 

r=0 

a{x) scalar, uniformly positive, and uniformly bounded. For the remainder 
of this section, we shall for notational convenience omit the subscript a, 
referring always to a-norms or -inner products unless otherwise specified. 
For later reference, we note the commutator relation 

(5.30) {f,g,) = -{f, + {aJa)f,g), 
and the related identities 

(5.31) (/,5/,) = -(l/2)(/, (5, + iaJa)S)f), 



(5.32) (/, (5/),) = (l/2)(/, (5, - (a./a)5)/), 

valid for symmetric operators S. 

By (Hl)(ii), we have that ^11(^51)"^ ^^^^ spectrum of uniform sign, 
without loss of generality negative, so that the similar matrix 

(^?,)-V2^n(A?,)-V2 = (A?,)-V2^^^(^o^)-i(^o^)i/2 
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has real, negative spectrum as well. (Recall, A^^ is symmetric negative 
definite as a principal minor of the symmetric negative definite matrix A^.) 
It follows that All itself is uniformly symmetric negative definite, i.e., 

(5.33) All <-e<o. 

Defining a, following Goodman |Go| . by the ODE 

(5.34) a^ = C,\W^\a, a(0) = 1, 

where C* > is a large constant to be chosen later, we have by ()5.33|1 

(5.35) ia^/a)Aii < -C,9\W^\. 

Note, because \ W,,\ < Ce"^!^!, that a is indeed positive and bounded from 
both zero and infinity, as the solution of the simple scalar exponential growth 
equation 1)5. 34(1 . 

Energy estimates for W: 

A^Wt + Aw, - (BW,), 

= -{A- A)W, + {{B - B)W,), 

-A—^5+{B 



(5.36) d5 \ 85 

- {A- A)—fb ^ [[B - B)—^] b 



where 
(5.37) 



db V ' 06 
A'^6-A'ipt-A^, + {Bip, 



A°:=A\W), A:=A{W), B := B{Wy, 
A°:=A^{W), A:=A{W), B := B{W); 



X 



(notice that we dropped tilde signs from A^,A and B in 1)5. 9|) and, with a 
slight abuse of notation, used them here differently). We want to write the 
right hand side of (|5.36|) in the form Aii + {Ai2)x + where M.i and A^2 
depending on W and "behaving well enough", and ^ is a remainder decaying 
fast enough. 

Beginning from the last line in (|5.36|) . we write 

~n dW^ ■ ~n n - dW^ dW^ ■ n - dW^ dW^ ■ 



47 



Now A%W+ip+^)^5 goes into ^ and (A^ -A°{W+^+^))^ 
goes into Mi. Notice that 



A°-A\W + + —-S)= / dA^W + ^ + —-5 + eW)d9W 
do Jq 86 

In a similar fashion AS'ipt and A(px each gives rise to a term which goes into 
M.I and the other which goes into ^. Similarly {Bipx)x comprises two terms, 
one of which is absorbed hj M2 and the other by ^. 
Using 

-dW^ - -dW^ -dW^ - -dW^ 
(5-38) dA^W. + A^ = {dB^Wx). + (S^). 

we can write the second and third lines of (|5.36|) in the form: 
(5.39) 

- - ~ - - -dW^ -dW^ 

{A - A)Wx-{{B - B)Wx)x + A— ^5 - {B—f)J 

00 00 

= (A-AiW + + —-6))W, 
00 

+ [a{W + ^ + -^b) -A- dA{^ + ^<5) j Wx + dA(pWx 

B - B{W + + -^S))Wx 

+ Ub{W + + -^b) -B- dBiip + -^6))WA 
-{dBipWx)^ 

the term in the second line of ()5.39|) goes into Mi, the third line and the 
fifth lines go into ^, and the fourth line goes into A^2- The fourth line of 
(|5.36j) can be dealt with in a similar way. 

To summarize, we were able to write equation (|5.36p in the form: 

(5.40) A'^Wt + AW, - (BW,), = Mi + (Mi), + t) 
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where M.i,M.2 are dependent on W; 

(5-41) _(^A-A{W + 0+ -Q^mx 

-^A-A{W + <f+—-5))W, 
oo 



-iA-AiW + ^+^5))^S. 



M2 ={B - B{W + (p+ -^S))ip^ 
(5-42) + - BiW + if + -^6))W, 

Using A = A + 0{(), Aj; = 0(|Wa;| + C) we can see: 

(5.43) \{dlW,dlMi)\ < C{dlW,W,dlW) + CC\W\lr + C\W\lr-, 

for r = 0, • • • ,3. M.2 has the block form: ^ J , hence using Young's 

yu M22J 

inequahty and the block structure of M2, 

(5.44) \{dlW,dl{M2h)\ < Ci,-'\w''\%r+tx\w''\lr+i 
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(5.45) 



for ^ arbitrarily small. is independent of W, 

- A^{W + <f + ^5)^t 

- A{W + ip + -^5)if^ 
+ {B{W + ^ + ^^5)^^)^ 

A{W + ^ + -^5) -A- dA{^ + + dA^W^ 

( - dw^ - - dW^ 

- UBiy/ + ^ + -^b) -B- dB{^ + -g^s))w, 

+ {dB^W,), 

-{AiW + ^ + ^5)-A)^5 
+ {{B{W + ^ + ^5)-A)^),6. 

^ has the good property that, with differentiation with respect to x, its rate 
of decay remains the same. Hence we transfer all the derivatives to ^: 

(5.46) \{DlW,Dli)\<C{\W\l, + \i\l,r) 

Notice that under the assumptions and definitions of Theorem 15. 5| 

(5.47) \^^.^s)\Hr <C[E^ + Q{t)){l + s)-l 

for any r and < s < t. Now the following lemma provides us with necessary 
bounds we need for |VF|//3: 

Lemma 5.7. Under the hypotheses of Theorem \5.5\ let Wq S H^, and sup- 
pose that, /or <t<T, both the supremum of\6\, \6\ and the W'^'°^ norm of 
the solution W = {w^,w^^y remain bounded by a sufficiently small constant 
C > 0. Then, for allO<t<T, 

(5.48) \W{t)\ls<C\Wmlse-'' + C [\-'-^'-^H\W\l, + \^\1,){t) dr. 

Jo 
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We first carry out a complete proof in the more straiglitforward case 
that the equations may be globally symmetrized , i.e., with conditions (Al)- 
(A3) replaced by the following global versions, indicating afterward by a few 
remarks the changes needed to carry out the proof in the general case. 
(Al') A^, Ai'' := i^, i° are symmetric, i° > 0. 

(A2') No eigenvector of J^^jdF^U) lies in the kernel of J2CjCkB^''iU), 

for all nonzero ^ G M'^. (Equivalently, no eigenvector of i-^^) ^(y^) 

lies in the kernel oi Y.CjCkB^^{W).) 

(A3') 



(5.49) 





6^'^ 



with ReY,(j^kV^{W) > 6'|Cp for some 6 > 0, for all W and all ^ € M.'^, and 
G = 0. 

To prove (|5.48|) . we carry out a series of successively higher order energy 

estimates of the type formalized by Kawashima [Kaw| and used extensively 
by K. Zumbrun et al (see [O] . also [ MaOl IMaZ.4p The origin of this 
approach goes back to jKan^ IMNi| in the context of gas dynamics; see, 
e.g., | HoZ.l| for further discussion/references. 

Let K denote the skew-symmetric matrix described in Lemma 15.111 as- 
sociated with A^, A, B, satisfying 



0\-l 



A + B>0. 



k{A 

Then, regarding A^ , K, we have 
(5.50) 

Al = dA'^{W)W^, = dK{W)W,, A, = dA{W)W^, 

io = dA\W)Wt, kt = dK{W)Wu At = dA{W)Wt, B, 



: dB{W)W^, 
dB{W)Wt. 



Now: 

dW^ - dW^ 

(5.51) \W^\ = \W, + W^ + ip:, + —^5\ <\W^\ + \W^\ + \^,\ + \—^5\ 

oo oo 

and, from the equation for A, 

(5.52) \Wt\<Ci\W,\ + \wii\). 
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Thus, in particular it follows that 



(5.53) 

1*^1) l^-rl' I^Sxl' l^j;;^;;!, \At\, \Bt\ 

<c{c + \u,\). 

In what follows, we shall need to keep careful track of the distinguished 
constant C^:. 
Computing 

(5.54) -{W,AW^) = {1/2){W, (i, + {a,/a)A)W) 

and expanding A = A + 0{C), A.j, = 0(1^^.1 + C), we obtain by 
the key property 
(5.55) 

-{W,AW,) = {l/2){w',{a^/a)Auw') 

+ 0(^{\W,\\WI\W\) + {ia,/a)\WlC\W\ + {w^'l)) 

< -{ae/3){\w^\,\w^\\w^\) + ccW\^ + cia)\w^^\^, 

by which we shall control transverse modes, provided C* is chosen suffi- 
ciently large, or, more generally, 
(5.56) 

-{d',w,Ad',w,) < -(a0/3)(|9,V|,|tr,pVl)+c^C|9Vl'+c^(c*)|9V 

Here and below, C(C=k) denotes a suitably large constant depending on 
C=K, while C denotes a fixed constant independent of C*: likewise, 0{-) 
indicates a bound independent of C*. 

Zeroth order "Friedrichs-type" estimate. We first perform a stan- 
dard, zeroth- and first-order "Friedrichs-type" estimate for symmetrizable 
hyperbolic systems [PriJ . Taking the a-inner product of W against (|5.4()|) , 
we obtain after rearrangement, integration by parts using (|5.3U|) - H5.31|1 . and 
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several applications of Young's inequality, the energy estimate 
(5.57) 

= -{W,AW,) + {W, (BW^U) + {W,Mi) + {W, {M2h) 
+ {W,0 + l{W,A^,W) 

= ^{W, (i, + {aJa)A)W) - {W^ - {a,/a)W,BW^) 
+ {W,Mi) - {W,M2) + {W,0 + l{W,A^tW) 

<-{W,,BW^) + C{a) J a(^{\W^,\ + \W,\)\W\^ 

+ \wi^\\W\i\W,\ + \W,\) + C\W\y + iiW^W^ + 

<-QW%^^c(c:)[\w\i,-,\i\y). 

Here, we used boundedness of and \Wx\ and also the inequalities 

(tCTTl) . (tCT^ and (tCTHl) . 

First order "Priedrichs-type" estimate. For first and higher deriva- 
tive estimates, it is crucial to make use of the favorable terms 1)5. 56() afforded 
by the introduction of a-weighted norms. Differentiating (|5.4nj) with respect 
to X, taking the a-inner product of Wx against the resulting equation, and 
substituting the result into the first term on the righthand side of 

(5.58) \{Wx,A^Wx)t = {Wx,{A^Wt)x) - {Wx,AlWt) + ^(VF,,i?W,), 
we obtain after various simplifications and integrations by parts: 

^{Wx,A'^Wx)t = -{Wx, iAWx)x) + {Wx, {BWx)xx) + {Wx, iMi)x) 

+ {Wx,{M2)x.) + {Wx,^x) 

-{Wx,AlWt) + ^{Wx,A^^Wx) 

=-{Wx,AWxx)-{Wx,AxWx) 

- {Wxx + {ax/a)Wx,BWxx + BxWx) 

+ {Wx, {Ml)x) - {Wx + {M2)x.) + {Wx,(x) 

- {Wx,AlWt) + ^{Wx,A^tWx). 

Estimating the first term on the righthand side of 1)5. 59() using (|5.56l) . k = 1, 
and substituting {A^)~^ times (|5.4()j) into the second to last term on the 
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righthand side of H5.59() . we obtain by (|5.53|) plus various applications of 

Young's inequality the next-order energy estimate: 

(5.60) 

+ CXCMW^'l + Cl^xl, (1^1 + \W,\)\W,\ + \wii\) 
+ C{i\W^\ + \wii\),\W^\{\W\ + \W^\)) 
+ C{W,,W^W,)+CC\W,\l2 



< -{e/2)\wii\l, - ic.e/i){m\wi\Awi\) + c{a.)cU\l, 
+ c{a)\w'j\i, + c{a){\w\L2 + \c\j,,), 

provided C=k is sufficiently large and C, fj, sufficiently small. 

First order "Kawashima-type" estimate. Next, we perform a "Kawashima- 
type" derivative estimate. Taking the o-inner product of against K(A^)~^ 
times (|5.4()j) . and noting that (integrating by parts, and using skew-symmetry 
of K) 
(5.61) 

= ^{W,,KWt)-^{Wt,KW,) 

-^{Wt, {K, + {aJa))W) + ]^{W.,,KtW) 
= {W^,KWt) + ^{W, {k, + {a^/a))Wt) + ^{W,,ktW), 
we obtain by calculations similar to the above the auxiliary energy estimate: 



-{W,,kW)t < -{W^,k{A'>)-^AW,) 

+ c{c.)\wi'f + cam + c + c)\wi\,\wi\) 

+ CC'\wii\' + C{CMW\L^ + \C\l.), 

where C > is an arbitrary constant arising through Young's inequal- 
ity. (Here, we have estimated term {AUx, {ax/ct)U) arising in the middle 
term of the righthand side of (t^lHl using (tOTl) by C(C*) J IW^WUI"^ < 

c{a)\u\i^.) 
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Combined, weighted H estimate. Choosing « ^ « 1, adding 
(|5.62() to the sum of H5.6U() times a suitably large positive constant M(C*, C) >> 

and ()5.57|) times M(C=k,C)^ ™d recalling 15.111 we obtain, finally, the 
combined first-order estimate 

(5 63) l{M{C,,Cf{W,A''W) + {W^,KW)+M{C,,0{W,,A''W^))^ 

< -oimh + Uiih) + c{c,) {\w\i, + , 

9 > 0, for any C(C)C'*) sufficiently small, and C*, C(C*) sufficiently 
large. 

Higher order estimates. Performing the same procedure on the twice- 
and thrice-differentiated versions of equation (|5.4()|) . we obtain, likewise, 
Friedrichs estimates 
(5.64) 

^{d^,wA'd^,w)t < -{e/2)\d^+W^\l, - icj/i){m\d>'mw'\) 

+ C7(c,)(ciavii. + i^vii^ + mHg'^ + \w\i, + ici^.,), 

and Kawashima estimates 
(5.65) 

^{d^,W,Kdr'W)t < -{d^,W,K{A'>)-^Ad^,W) 

+ c{a)\diw''\i, + cam + c + cmw'i 

+ CC'm^Wll, + C7(C7.)(|T4^.|^.-. + \W\l, + 

for g = 2, 3, provided C(C, C=k) are sufficiently small, and C=k, C(C=k) are 
sufficiently large. The calculations are similar to those carried out already; 
see also the closely related calculations of Appendix A, |MaZ.2j . 

Final estimate. Adding M{C^,Cf times ((EMI, M{C^X) times 
and ()5.65() . with q = 2, where M is chosen still larger if necessary, we obtain 

(5.66) 

^ (m(C,, C)^(Ty, i^WO + M{C,,Cf{W,,KW) + M{a,cf{w.,A''w,) 
+ {dlW, Kd^W) + M{a,C){dlW, A'dlW))^ 

< -eimir + \wi'\i,) + c{c,) [\w\i,) + iei^4) . 

Adding now M{C^,C,f times (jSSZI), M(C=^,C) times (jSSU, and ((EM]), 
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with g = 3, we obtain the final higher-order estimate 
(5.67) 

+ M{c,,cf{dlw,Kd^w) + M{c,,cf{dlw,A^dlw) 

+ {dlW,KdlW) +M{C^,C){dlW,A^dlW)^^ 

< -oimi^ + \wi'\i,) + c{c,) [\w\i, + lei^e) . 

<-e\w\]j,+c{c,){\w\l,) + \i\]j,). 

Denoting 
(5.68) 

£{W) := ^(^M{C.,Cf{W,A''W) + M{C.,C)''{W^,KW) 
+M(a, Cf{W^, A'^W^) + M{a,Cf{dlW, Kd^W) + 
MiCCfidlW, A'dlW) + {d^W, kdlW) + M(C,X){dlW, A'^dlW)) , 

we have by Young's inequahty that £^^'^ is equivalent to norms and 
H^, hence (|5.6.S|1 yields 

£,<-0^£ + c{a){\w\l,) + \c\jje), 

from which we conclude, 

£{t) < e~'^'£{0) + C{a) f e"^2(*-^) {\W\l,) + l^l^e) {s) ds. 

Jo 

This is equivalent to (|5.48jl . 

The general case. It remains only to discuss the general case that 
hypotheses (A1)-(A3) hold as stated and not everywhere along the profile, 
with G possibly nonzero. These generalizations requires only a few simple 
observations. The first is that we may express matrix A in (|5.4fljl as 

i=A+(i«'.i+o(o(i) 

where A is a symmetric matrix obeying the same derivative bounds as de- 
scribed for A, identical to A in the 11 block and obtained in other blocks 
jk by smoothly interpolating over a bounded interval [—R, +R] between 
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A{W-oo)jk and A{W+oo)jk- Replacing A by A in the qth order Friedrichs- 
type bounds above, we find that the resulting error terms may be expressed 
as (integrating by parts if necessary) 

{d^,o{m+c)\w\,\d^+w^\) 

plus lower-order terms, hence absorbed using Young's inequality to recover 
the same Friedrichs-type estimates obtained in the previous case. Thus, we 
may relax (Al') to (Al). 

The second observation is that, because of the favorable terms 

-{ae/A){m\d^,w'\,\diw'\) 

occurring in the righthand sides of the Friedrichs-type estimates, we need 
the Kawashima-type bound only to control the contribution to IdxW^l"^ com- 
ing from X near ±00; more precisely, we require from this estimate only a 
favorable term 

-e{{i-o{\Wx\ + C + C))\d>'\Mw'\) 

rather than —OlOxW^l"^ as in ()5.62|) and (|5.65l) . But, this may easily be 
obtained by substituting for K a skew-symmetric matrix- valued function K 
defined to be identically equal to -^(-|-oo) and ^(—00) for |2;| > R, and 
smoothly interpolating between ^(±00) on [—R,+R], and using the fact 
that 

{K{A^)-^A + B)^ >9>0, 

hence 

{K{A^)-'A + B)> 6(1 - 0{\Wx\ + 0). 

Thus, we may relax (A2') to (A2). 

Finally, notice that the term G — G in the perturbation equation may 
be Taylor expanded as 

{~g{Wx,Ux) + g{Wx,Wx\) ^ {o{\Wx\^)) 

The first, linear term on the righthand side may be grouped with term 
A^Wx and treated in the same way, since it decays at plus and minus spa- 
tial infinity and vanishes in the 1-1 block. The (0,O(|Wj;P) nonlinear term 
may be treated as other source terms in the energy estimates Specifically, 
the worst-case terms {d^W, KdlO{\Wx\'^)) and {d^W,d^{0,O{\Wx\'^))) = 
{d^w^^ ,d'^0{\Wx\'^)) may be bounded, respectively, by |VF|vi/2,oo|H^|^3 and 
|VK|iy2,oo|w-'^-'^|ji^4|I^|ji^3. Thus, we may relax (A3') to (A3), completing the 
proof of the general case (Al)-(A3) and the theorem. □ 
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